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INFLUENCE OF PARTICLE DRAG COEFFICIENT ON 
PARTICLE MOTION IN HIGH-SPEED FLOW WITH TYPICAL 
LASER VELOCIMETER APPLICATIONS 

Michael J. Walsh 
Langley Research Center 

SUMMARY 

The effect of using different particle drag coefficient Cq equations for comput- 
ing the velocity of seeded particles in high-speed gas flows has been investigated. The 
Cj) equations investigated include the Stokes equation, a second incompressible equation 
valid for higher relative Reynolds numbers, and six more detailed equations that account 
for the effects of compressibility together with the effects of relative Reynolds numbers 
greater than one. The flows investigated are center-line nozzle flows, normal shocks, 
and oblique shocks for free-stream Mach numbers of 1.6 to 6 and stagnation pressures 
of 1 and 3.4 atmospheres. Particle sizes range from 0.5 to 10 jum. These flows were 
selected because of their similarity to flows encountered in previous laser velocimeter 
studies in supersonic flows. The accuracy of the data on which the empirical Cq equa- 
tions are based was also investigated. The net result is an empirical Cq equation 
based on the latest sphere Cd data for the low relative Mach number and Reynolds 
number conditions that are encountered in supersonic flows. This new Cq equation 
is also used to examine the effect of gas density on the relaxation length behind a Mach 6 
normal shock for stagnation pressures from 0.3446 to 20.7 MN/m^ (3.4 to 204 atm). 

(1 atmosphere equals 101.3 kN/m^.) 


INTRODUCTION 

The laser velocimeter (LV) is an instrument used to determine flow-field velocities 
by measuring the velocities of particles seeded in high-speed gas flows. Problems occur 
in the application of the LV to flows where large velocity gradients are present. The 
presence of such gradients in many supersonic flows creates a situation where the parti- 
cles cannot accelerate or decelerate as rapidly as the gas does. The resulting difference 
between the gas velocity and particle velocity is called the particle "velocity lag." When 
using an LV system in a supersonic flow, the investigator must consider this particle 
velocity-lag error. 



To measure gas velocity in a supersonic flow, limitation of particle diameters to 
approximately 1 jim may be necessary to minimize particle velocity lag. However, the 
Mie scattering criteria (ref. 1) and the characteristics of the LV system (as discussed in 
ref. 2) may limit the minimum particle size used in an LV application. Thus, particle- 
motion studies are necessary to determine the size of particles required to follow the gas 
flow; when the Mie scattering criteria limits the minimum size, these calculations are 
necessary to determine the particle velocity lag. The particle motion is governed by the 
flow-field properties, the particle properties, and the particle drag coefficient Cj). A 
search of the literature in the area of high-speed flows reveals that a number of equations 
have been used for particle-motion studies, ranging from the simplest Stokes Cp) equa- 
tion (ref. 3) to the more detailed equations of references 4 to 10. To date only two studies 
have examined the influence of Cj) on particle-motion calculations. Maxwell and 
Seasholtz (ref. 11) used two incompressible Cp) equations and the equation of Carlson 
and Hoglund (ref. 6) to compare the particle -motion calculations behind Mach 1 to Mach 2 
normal shocks. In reference 10, Korkan, Petrie, and Bodonyi used the equations of ref- 
erences 4, 6, and 7 to compare particle -motion calculations for uniform flow, Prandtl- 
Meyer expansions, and oblique shocks in a Mach 5 flow. 

The objectives of this investigation are: (1) to determine how well the Cp equa- 

tions of references 4 to 10 predict the latest available sphere Cp) data and, if necessary, 
to improve the accuracy of sphere Cp) "predictions"; (2) to determine the influence of 
Cp) on particle-motion calculations for various high-speed gas flows using the Cp> equa- 
tions of references 4 to 10; and (3) to examine the influence of gas density on particle 
motion behind a normal shock. 


SYMBOLS 


Cp) drag coefficient 

D diameter, pm 

M Mach number 

^Kn Knudsen number based on particle diameter 

NRe,r relative Reynolds number 

p pressure, N/m^ 

2 
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R specific gas constant, J/kg-K 

T temperature, K 

t time, s 

V velocity, m/s 

X distance, m 

y ratio of specific heats for gas 

9 percent velocity lag (eq. (14)) 

X relaxation length, distance from shock to point v^here 0 = 1 percent, cm 

jU viscosity, N-s/m^ 

p density, 919.5 kg/m^ 

Subscripts: 


c 

continuum 

comp 

compressible 

FM 

free molecular 

g 

gas 

inc 

incompressible 

max 

maximum 

P 

particle 

r 

relative 

X 

X -component 


3 


y 


y -component 


0 tunnel stagnation conditions 

oo free stream 


An arrow over a symbol represents a vector quantity, 

DISCUSSION 


Governing Equation for Particle Motion 

When applying laser velocimeter (LV) systems to gas flows where particle velocity 
lag is significant, the motion of the particle becomes important. In these LV applications 
the particle mass densities are typically much greater than that of the gas. The govern- 
ing equation for a spherical particle traveling in a fluid where the particle mass density 
is much greater than the gas density is given by Soo (ref. 3) as 


dVp 

dt 



( 1 ) 


By defining a relative Reynolds number based on the velocity difference between the gas 
and the particle as 


N- 


Re,r 



( 2 ) 


equation (1) becomes 


^-3 

dt 4 


^D^Re,r^g 



(3) 


Thus, the particle acceleration depends on Nj^g p, flow -field properties, particle mass 
density, particle size, velocity lag, and Cj). 

As discussed in references 4 and 5, equation (3) can be treated in component form 
for two-dimensional flows as 


4 



( 4 ) 


3 ^E)N^Re,rMg(Vg,x " ^p,x) 


dt 




and 


dVp,y _ 3 <^D^Re,rMg(Vg,y - Vp,y) 
^ PpDp^ 


(5) 


To complete the system of governing differential equations, two additional differential 
equations are obtained from the definition of velocity: 


dX = V 


dy TT 

^P"Vp,y 


(6) 

( 7 ) 


Thus, a numerical solution of equations (4) to (7) gives the particle motion through any 
known two-dimensional flow field, provided the components of initial particle velocity and 
position, particle properties, and gas properties are known. Various mathematic repre- 
sentations of the Cd term appearing in equations (4) to (5) are discussed in the follow- 
ing section. 


Review of Methods for Calculating Cq 

A search of the literature concerning high-speed flows revealed that theoretical 
solutions for Cp exist only for low relative Reynolds number (Nrq r < 1), incompressi- 
ble flows, and free -molecular flows. For these Nrq ^ values the Stokes Cp equation 
from reference 3 


Cd = 


24 

^Re,r 


( 8 ) 


is valid. For free-molecular flows assuming diffuse reflection, Emmons (ref. 12) gives 
the following Cd equation: 


Cd,fm = (l + 2s2) exp erf (s) + 2 V? 


1/2 


(9a) 


5 



where 


s = 



(9b) 


and where the relative Mach number is defined as 





(9c) 


For incompressible flows where Nrq r ^ there is a significant amount of experi- 
mental data. Empirical equations have been derived which predict these data for limited 
ranges of NRg,r- The empirical Cq equation given by Torobin and Gauvin (ref. 13) is 


Cd = 



( 10 ) 


Equation (10) gives good predictions of the incompressible steady-state sphere Cp data 
tabulated by Perry (ref. 14) for < 200. 

For many LV applications in high-speed flows, the dependence of the drag coefficient 
on Mj. (compressibility effect) as well as > 1 must be accounted for. The avail- 

able Cp) equations that are applicable to high-speed flow are the ones used by Cuddihy 
et al. (refs. 4 and 5), Carlson and Hoglund (ref. 6), Crowe (ref. 7), Crowe et al. (ref. 8), 
Waldman (ref. 9), and Korkan et al. (ref. 10). These equations are longer and more 
detailed than equations (8) to (10) given above, and are given in the appendix. 

To date LV researchers have used several different Cj) equations for particle - 
motion calculations. Asher et al. (ref. 15), assuming small velocity lags, use Stokes drag 
equation (8) to determine turbulence velocity spectra. In their early work at Mach 3, 

Yanta et al. (ref. 16) used the following incompressible Cj) equation from reference 17: 


Cd = 28N^°-S/ + 0.48 (11) 

In later work at the same Mach number, Yanta (refs. 18 and 19) uses the Cd equation 
from references 4 and 5. In an LV application study at Mach 5, Meyers and Walsh (ref. 2) 
also use the method from references 4 and 5. Maxwell and Seasholtz (ref. 11) examined 
particle motion through a normal shock for a free -stream Mach number 1.6, and Maxwell 
(ref. 20) computed particle motion through turbomachinery using the Cd of Carlson and 


6 



Hoglund, Morse et al. (ref. 21) examined particle motion in a Mach 5 flow using an 
earlier version of Crowe’s method. 

In considering the various Cp) equations that have been used in past LV studies, it 
is important to determine which equation is the most applicable to LV studies. The next 
sections examine the available experimental sphere Cj^ data and compare the predictions 
of the Cd equations in references 4 to 10 to the sphere Cp) data to determine the most 
suitable equation for LV use. 

Experimental Data 

A review of the literature which studies drag coefficient of spheres indicates that 
Bailey and Hiatt (ref. 22) provide the most extensive sphere Cj) data available. The 
Mr < 2 and NRe,r 200 range covered by the experimental Cp) data of Bailey and 
Hiatt (ref. 22) and of references 23 to 27 are shown in figure 1. Figure 1 indicates that 
the Bailey and Hiatt data cover a significant part of the Nfje r 200 range and the 
Mr < 2 range that was not covered by previous data. Figure 2 shows that the variation 
of the previous drag coefficient data of references 23 to 29 and the Bailey and Hiatt data 
is less than 2 percent in the NRg r I'ange 10^ to 10® and as much as 13 percent at the 
lower Nrq r (NRg r = 30, Mr = 2). 

Bailey (ref. 30) has examined some experimental factors, such as turbulence and 
model support interference, that affect the sphere Cjy measurements; his study has 
shown that most of the early sphere drag data are in reasonable agreement with the Bailey 
and Hiatt data if these factors are accounted for. Zarin has observed that turbulence has 
little effect on Cp) for < 100 if the turbulence intensities are below 3 percent. 

Aroesty (ref. 31) and Sherman (ref. 32) have also noted that model support interference 
may cause large errors in Cj) measurements. The only data shown in figures 1 and 2 
that are free of model support interference and turbulence effects are the data of Bailey 
and Hiatt, of Zarin, and of Goin and Lawrence. The data of Goin and Lawrence are limited 
but show excellent agreement with the data of Bailey and Hiatt. As shown in figure 1, the 
data of Zarin cover a lower Mp range than the data of Bailey and Hiatt. Therefore, the 
combined data of Zarin and of Bailey and Hiatt provide the most complete and accurate 
coverage of the Mp and p range encountered in partical -motion studies connected 

with LV systems. The following section compares the Cp) equations of references 4 
to 10 to this sphere Cj) data. 

Selection of the Most Applicable Cp> for LV Studies 

The equations of references 9 and 10 are the only ones of references 4 to 10 that 
were published after the Bailey and Hiatt sphere Cp) data became available. Without 
referencing any new data, Waldman (ref. 9) modified the earlier method of Crowe (ref. 7); 
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Korkan, Petrie, and Bodonyi (ref. 10) used the data of Bailey and Hiatt (ref. 22). A com- 
parison of the equations used in references 4 and 5 with the equations in reference 10 (see 
appendix), together with a comparison of the parameter values in tables I and II, indicate 
that there is little difference between the two methods except for the temperature correc- 
tion used in reference 10 and differences in the continuum Cp values Cp q . The value 
of Cj 5 0 is the value of Cq at large Nj^e,r- Therefore, the method used in refer- 
ence 10 has only modified the high NRe,r predictions of references 4 and 5. 

Figure 3 compares the predictions of the six Cq equations of references 4 to 10 
with experimental sphere Cjj data to determine the best Cj) equation available. For 
the comparisons in figure 3, the particle or sphere temperature is assumed to be equal 
to the gas temperature (the condition of the experimental data). As noted earlier, refer- 
ence 10 used only the higher Nj^e,r reference 22. Figure 3 shows that the method 

of Korkan et al. fails to predict the data of reference 22 for p < 100. Figure 3 also 

shows that the other methods (refs. 4 to 9) do not give good predictions of low Npg p and 
Mr sphere Cd data. 


Better Predictions for Cp 

In this section the method used by Cuddihy, Beckwith, and Schroeder (refs. 4 and 5) 
is modified to give better predictions of the experimental sphere Cjj data of refer- 
ences 22 and 26. For Mp ^ 0.5, the method given in references 4 and 5 used the fol- 
lowing equation: 


= Cd,C + (Cd,FM - Cd,c) exp 


-A(NRe,r)^ 


( 12 ) 


where Cp)^c are the continuum and the free -molecular values of Cp), 

respectively. The parameters A and N are functions of Mp and are selected to fit 
experimental sphere Cp) data. With suitable mathematical operations, equation (12) 
becomes 


loge 


loge pD,FM - Cd,c) - l°ge pD - Cd,c) = loge A + N logg N^g^p 


(13) 


which is the equation of a straight line in the coordinates 




Cd,fm - Cd,c\ 
Cd - Cd,C / 


; logg NRg^ 


r 
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A least squares fit of available experimental sphere Cj) data then yields the value of A 
and N, 

This paper uses equation (12) for Mj- ^0.1 and adjusts the parameters A and 
N to fit the experimental Cy) data of references 22 and 26. Also, the Cp> values are 
required to approach the incompressible Cd values of equation (10) as Mj. approaches 
0.1 and compressibility becomes negligible. The new values for the parameters Cj) q, 
Cj3, A, and N in equation (12) are given in table HI. The values for Cp) q are 
obtained from the high data of reference 22. These values were selected on the 

basis that Cj) q is the value Cj) approaches at large Thus, Cq c can be 

tabulated as a function of Mr- For LV applications, the gas molecules would reflect from 
the particles diffusively; therefore, the Cj) values are determined by use of equa- 
tion (9) which assumes diffuse reflection. The values of used by references 4 

and 5 for their particle -motion studies were the values of Cd,FM, assuming specular 
reflection. 

Figure 4 indicates that the method presented here gives excellent predictions of the 
experimental Cp> data. These equations do not account for differences between particle 
temperature and gas temperature since there is no reliable experimental data on which to 
base a temperature correction for relative Mach numbers less than 2, 

The small region of extrapolation in this approach (Mr <1.0; NR0 ^r ^ 40) can be 
further reduced by introducing into the experimental data base any new low relative Mach 
number and Reynolds number data that become available. Based on comparisons, it has 
been found that equation (12), with the parameter values given in table III, gives the best 
predictions of sphere Cp) data; therefore, it is more applicable to particle -motion stud- 
ies in connection with LV systems. The next section examines the influence of Cq on 
particle motion for a limited number of high-speed gas flows. 

Effect of Cj) on Particle -Motion Calculations 

As mentioned earlier, particle motion is sensitive to the size, mass density, and 
initial velocity of the particle as well as to the gas flow -field properties. This report 
is concerned only with the influence of Cjy on the particle-motion calculations. This 
influence is determined by examining a limited number of high-speed flows that may be 
encountered in LV applications. 

Table IV lists the stagnation pressure, stagnation temperature, gas velocities, max- 
imum Mp and encountered by the particles during the calculation, and the loca- 

tion in the flow where the particle -velocity calculations are compared. Table V lists the 
dimensions of the nozzles that are important to the particle -motion calculations. The 
particle-mass density for all the test cases was equal to 919.5 kg/m^. It should be noted 
that the location of comparison is arbitrary and that the particle velocity lag varies with 
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distance. Thus, the differences between the particle velocities predicted using the various 
Cd equations may increase or decrease with a change in the comparison location. For 
all the test cases where the particles traversed a shock, the initial particle velocity was 
assumed to be equal to the free -stream gas velocity ahead of the shock. For the Mach 5 
and Mach 6 center-line test cases, the particles were injected upstream of the nozzle 
throat at a velocity equal to the gas velocity at that point, 30,48 m/s. For the Mach 3 test 
case, the particles were injected at the throat of the nozzle at the sonic gas velocity. As 
shown in tables IV and V, the comparison locations for the particles flowing along the cen- 
ter line of the nozzle were selected so as to be on or downstream of the nozzle exit. For 
the test cases involving particles traversing a shock (see figs. 5 to 7), a comparison loca- 
tion was arbitrarily chosen to be 1,27 cm behind the shock. For the oblique shock, the 
two-dimensional particle motion was calculated, whereas one-dimensional particle motion 
was calculated for the normal shock and the center-line nozzle test cases. 

Figures 5 to 7 indicate the percent velocity lag calculated by using various Cp) 
equations, where the percent velocity lag of the x-component is defined as 


0x = Percent velocity lag of x-component = 100 


Vg,x - Vp,x 


V, 




and that of the y-component as. 


6y s Percent velocity lag of y-component = 100 


^g,y ■ ^p,y 


g,y 


(14a) 


(14b) 


Figures 5 and 6 examine the influence of Cp) on particle -motion calculations using both 
the methods from references 4 to 10 that consider compressibility and the present method; 
figure 7 examines the incompressible Cp) equations (8) and (10) in comparison to equa- 
tions used in the present method. For LV applications in high-speed flows, it is expected 
that particle sizes would be in the 0.5- to 2-|Um range. Figures 5 to 7 show that even in 
this particle-size range, the percent velocity lag predicted by the previous Cp) equa- 
tions may vary considerably for certain flows. For example, a 0,5-/im particle passing 
through a normal shock in a Mach 3 flow has a particle velocity lag somewhere between 
0 and 80 percent, 1.27 cm behind the shock, depending on the Cp) method used. This 
large uncertainty in the velocity lag is caused by variations in the Cp) predictions. The 
present method gives better predictions of the Cp> sphere drag data and, therefore, bet- 
ter predictions of the particle motion. Figures 5 to 7 indicate the importance of using an 
accurate Cp> equation in the regions behind normal shocks and in the y-component of 
particle-velocity calculations for 5° and 10° oblique shocks in the Mach 3 to Mach 6 flows 
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for the stagnation conditions listed in table IV. As mentioned earlier, the flow fields exam- 
ined in this paper are limited in number and were selected only to determine whether the 
Cj) equation used could have an important influence on particle-motion calculations. It 
must be noted that there may be other flow-field conditions or measuring locations where 
inaccurate Cd methods may lead to larger errors in particle- motion predictions than 
determined in this report. 

Figures 5 to 7 show that large percent velocity lags occur when particles pass 
through a normal shock. The previous discussion has been limited to examining particle- 
velocity predictions behind normal shocks at 1.27 cm behind the shock. In supersonic 
flows where normal shocks may occur, it is important to determine how close to the shock 
accurate LV measurements can be made. Stokes Cp) equation is often used in particle - 
motion calculations since it has an analytical solution for constant velocity fields. (See 
ref. 17.) Figures 8 to 11 present the percent velocity lag as a function of distance for the 
Mach 1.6 to Mach 6 normal shocks and stagnation conditions examined earlier. The cal- 
culations were performed by use of the present method and Stokes Cp> equation. The 
figures show that the differences in velocity-lag predictions do vary with the comparison 
location as mentioned earlier. 

Often, the distance behind the shock where particle velocity lag has decreased to a 
specified percent is needed. This distance is defined as the relaxation length (in cm). 

The specified percent used in this paper is 1 percent. Figure 12 gives the relaxation 
lengths for the Mach 1.6 to Mach 6 normal shocks and stagnation conditions examined pre- 
viously. There are considerable differences in the calculated relaxation lengths as the 
particle size increases. The Stokes method predicts greater relaxation lengths than the 
present one for Mach 1.6 and Mach 5 (except for Dp< 3.5 jum) normal shocks; however, 
the situation is reversed for the Mach 3 and Mach 6 normal shocks. This reversal in the 
trends of X with Dp is caused by the local density which is dependent on the stagnation 
temperature and pressure and the free -stream Mach number. 

The effect of stagnation pressure on the calculation of relaxation length behind a 
Mach 6 normal shock is shown in figure 13. At low stagnation pressures, Stokes Cp 
equation predicts smaller relaxation lengths than the method used here. As the stagnation 
pressure increases, the present method predicts smaller relaxation lengths; however, the 
Stokes calculations are not affected by changes in stagnation pressure. If Stokes Cq 
equation (8) is substituted into equation (4), the following equation is obtained for the par- 
ticle motion behind a normal shock: 


dVv 


dt 


18^^g(Vg,x - Vp,x) 

PpDp2 


(15) 
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Equation (15) clarifies the results shown In figures 12 and 13. Particle -motion calcula- 
tions using Stokes Cj) equation are independent of changes in the gas density or pressure 
and are affected by temperature changes only through the gas viscosity jUg. Reference 33 
indicated that stagnation temperature and free-stream Mach number have a much smaller 
effect than stagnation pressure on the errors in using Stokes Cj^ equation for particle- 
motion calculations. In summary, Stokes equation may be a good simple Cj) equation 
for low stagnation pressure calculations, but it may lead to large errors at large stagna- 
tion pressures. 


CONCLUDING REMARKS 

A number of particle drag coefficient Cj) equations are available in the literature 
on sphere drag coefficients. Particle motion in a number of supersonic flows was exam- 
ined using these Cp) equations. There was little difference in the velocity predictions 
for center-line nozzle flows and the horizontal components of particle velocities behind 
the oblique shocks. However, large variations occurred in the velocity predictions for 
particles passing through normal shocks and in the vertical components of particle veloc- 
ities behind oblique shocks. 

Available sphere Cj) data provided an experimental data base that was used to 
evaluate the predictions of the various Cp) equations. It was found that none of the avail- 
able methods could accurately predict the low relative Mach number and relative Reynolds 
number (m^ and Nj^e,r) experimental sphere Cj) data. 

A modified version of the method used by Cuddihy, Beckwith, and Schroeder has been 
developed to give accurate predictions of the available low Nj^^^p and Mj> experimental 
sphere Cj) data. The new method was then used to demonstrate the importance of an 
accurate Cp) equation for particle -motion calculations behind normal shocks and calcu- 
lations of the y-component of particle velocity behind 5^ and 10^ oblique shocks in Mach 3 
to Mach 6 flows. Finally, it was determined that errors in using Stokes Cp) equation for 
the calculation of relaxation lengths behind normal shocks were extremely sensitive to the 
stagnation pressure of the free stream. These errors increased as the stagnation pres- 
sure increased. 

In conclusion, an accurate Cj) method is needed for particle- motion studies 
required in laser velocimeter diagnostics in high-speed flows since wide variations can 
occur in particle-velocity calculations for some supersonic flows if the Cj) method 
used does not accurately predict the sphere Cj) data. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
January 5, 1976 
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APPENDIX 


DETAILS OF PREVIOUS DRAG COEFFICIENT METHODS 

Any evalviation of the Cj) methods that account for compressibility requires an 
examination of the details of the various equations. This appendix discusses the details 
of the individual Cj) methods of references 4 to 10, and similarities are noted (to explain 
why all existing prediction methods fail to provide accurate predictions of low Mj. and 
NRe,r experimental sphere Cj) data). 

Method of Cuddihy, Beckwith, and Schroeder 

For Mj. greater than or equal to 0.5, the Cjy equation used by Cuddihy, Beckwith, 
and Schroeder (refs. 4 and 5) is 


Cd = ^D,C + (Cd,FM - ^D,c) exp 


-A(NRe,r)” 


(Al) 


The values for the continuum drag coefficient Cq^C for Mj. = 1.6 to 9.7 are taken from 
the ballistics range data of May and Witt (ref. 29) and Hodges (ref. 34) for Nj^e,r “ 
to 10®. The free-molecular drag coefficient values Cj) fm taken from Emmons 
(ref. 12). The data used for the evaluation of the A and N parameters are taken from 
the experimental data of Aroesty (ref. 23), Sreekanth (ref. 24), and Ashkenas (ref. 25). 

The values for the parameters Cj) c, Cj) pjyj. A, and N are listed in table I. 

For Mr < 0.5, the Cp equation used by Cuddihy, Beckwith, and Schroeder (refs. 4 
and 5) is based on the high Nj^g r data (10'^ to 10®) of Charters and Thomas (ref. 28) and 
the incompressible steady -state Cp curve given by Rouse (ref. 35). The resulting equa- 
tion is 


Cd 


^^51.1Mr 


N 


Re,r 


l + 0.256Mrfc + ^^-— r 


^Re,i 


(A2a) 


where 


C = — — — + 0.4 + 1.6 exp 
NRe,r 


-(o.028N^®|) 


(A2b) 
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APPENDIX 


It is noted that the data of Bailey and Hiatt and of Zarin was not available when this 
Cj) equation was developed. Figure 1 indicates the large amount of extrapolation neces- 
sary to predict Cj) at low values of NRg,r ^r equations (Al) and (A2) are used. 

Method of Carlson and Hoglund 

Carlson and Hoglund (ref. 6) approach the development of a Cd equation by modi- 
fying the Stokes Cj) equation (8) for NRe,r > 1 as done in reference 13: 

CD,inc = + 0.15N^e®f ) (A3) 

^^Re,r 

Next, equation (A3) is corrected for compressibility effects based on an empirical corre- 
lation of the high NRe r sphere Cp) data presented by Hoe rner (ref. 36): 



Finally, using the work of Millikan (ref. 37) and the experimental data of Stalder and 
Zurick (ref. 38), Carlson and Hoglund modify equation (A4) to account for rarefied flow 
effects: 



One disadvantage of equation (A5) is that the constants 3.82 and 1.28 are evaluated to give 
the correct rarefied flow Cp) for Mj. = 0.5. For any other Mj., equation (A5) does not 
give the correct free -molecular limit of Cp). 

Carlson and Hoglund could only compare the predictions of equation (A 5) with Mach 2 
experimental data since the lower Mj. and Nj^g j. data of Bailey and Hiatt and of Zarin 
were not available in 1964. 

Method of Crowe 

A later paper by Crowe (ref. 7) compared the predictions of the equation of Carlson 
and Hoglund with Mach 3 experimental data and found poor agreement. Crowe then devel- 
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oped a Cd method that predicts Mach 3 experimental data. The equation was formulated 
by Crowe to give the correct value of the Cp) at = 0 and M^. = 2. The resulting 
equation is 


Cd = (CD,inc - 2) exp 




hfMr) 


Npe,r\ 

2Mr / 


+ 2 


(A6a) . 


where 


log 10 g(NRe,r) 


1.25 


1 + 


tanh (0.77 logio NRg,r 



(A6b) 


/T \^/2 

h(Mr) = 2.3 - 1.7 — - 2.3 tanh (1.17 logjQ ^r) 

VsJ 


(A6c) 


Both the methods of Crowe and of Carlson and Hoglund approach the incompressible 
Cj) values and account for compressibility effects by correlating the compressible high 
NRe,r data of Hoerner. 

Method of Crowe, Babcock, and Willoughby 

In a later paper Crowe, Babcock, and Willoughby (ref. 8), noting that equation (A6) 
exhibited an unlikely inflection point, developed a new Cd equation that depended more 
on experimental sphere Cd data. The basic equation of Crowe, Babcock, and Willoughby 
is 


Cd = Cd,c + (Cd,fm - Cd,c)Cd 


A comparison of equations (Al) and (A7) shows that 


Cd = exp 


-A(NRe,r)^ 


(A7) 


(A8) 


The main difference between the two equations is that the method used by Cuddihy, 
Beckwith, and Schroeder tabulated the parameters Cd^C) Cd,FM> A, and N as given 
in table I, and Crowe, Babcock, and Willoughby provided equations for the various param- 
eters Cd C> Cd,FM> Cd- Also, Crowe, Babcock, and Willoughby were able to base 
their equation on additional experimental sphere Cd data that had become available since 
the development of the method used by Cuddihy, Beckwith, and Schroeder. The Cd,FM 
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values for equation (A7) were obtained from equation (9) given by Emmons, and the 
data were taken from an empirical correlation of the compressible Cj) data of Hoerher 
as given below 


Cd^C = 0-®® + 0.26 tanh (2 logg Mr) + 0.17 exp 


. 2 . 5 ( 12 ^' 


(A9) 


A comparison of equations (A5), (A6), and (A9), makes it difficult to see any similarities 
even though Hoerner's data were used for the correlations in all three equations. The 
equation for Cp in equation (A7) was determined by correlating the experimental data 
of Aroesty (ref. 23), May and Witt (ref. 29), Sims (ref. 39), Sivier and Nicholls (ref. 40) 
Zarin (ref. 26), and Millikan (ref. 37) among others. The resulting equation for Cq is 


Cd = 


K 

K + 1 



exp 


"®^Kn 



(A 10a) 


where 


K = N^ exp 


(AlOb) 


(^D, inc " 0.4)NRg^r 
8 


(AlOc) 


Method of Waldman 

Waldman (ref. 9) modified Crowe's method to account for the transonic drag rise at 
high Reynolds number without making any reference to additional data. The resulting 
equation used by Waldman is 


Cd = 


CD,inc - 2 + 0.463F(Mr) 


exp 


Mr 


NRi 


e,r 


= (NRe,r) 


+ R(Mr) exp f-0.5 + 2 


(Alla) 


where 


F(Mr) = exp 


0.863 

Mr^ 


1.163\ 
Mr^ j 


(Allb) 
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z(NRe,r) = 3.632 exp 


( ^ 



R(Mr) = 


1 

1.183Mr 


4.65 ^ 
Mr + 1 



(Allc) 


(A lid) 


Method of Korkan, Petrie, and Bodonyi 

The latest method found in the literature is the method of Korkan, Petrie, and 
Bodonyi (ref. 10) who compared the drag coefficient equations of Cuddihy, Beckwith, and 
Schroeder, of Carlson and Hoglund, and of Crowe to the 1971 sphere drag data of Bailey 
and Hiatt. They determined that the method used by Cuddihy, Beckwith, and Schroeder 
gives the best prediction of the Bailey and Hiatt data. The equations as presented in ref- 
erence 10 are slightly different from the equations used by Cuddihy, Beckwith, and 
Schroeder. Private communication with the author of reference 10 confirmed that for 
Mr <0.5 the drag coefficient equation should be identical to equation (A2) originally used 
by Cuddihy, Beckwith, and Schroeder. Also, for Mr = 0.5 the drag coefficient equation 
should be identical to equation (Al) used by Cuddihy, Beckwith, and Schroeder except for 
the addition of a temperature correction term, ACd(T), and the Cq^C parameters in 
equation (Al). The new C[> q values are listed in table II. Equation (Al) used in ref- 
erence 10 with a temperature correction term is given below: 


Cd = Cd,C + (^C,FM - ^D,c) exp -A(NRg,r)^ 


+ ACd(T) 


where for Mr > 1 


- 1 


ACd(T) 


1.183Mr 


0.142 + 12.22 - 0.597,|^) exp 

Mr 


(A 12a) 


(A 12b) 


and for Mr < 1 

ACd(T) = 0 (A12c) 

In summary, the only improvement of the method of reference 10 over the one used 
in references 4 and 5 would be the predicted Cj) values for high values of Npjg j.. 
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TABLE I,- PARAMETERS IN DRAG COEFFICIENT EXPRESSION 

[From refs. 4 and 5] 


Mr 

Cd,c 

Cd,FM 

A 

N 

0.5 

0.520 

7.80 

0.315 

0.410 

.6 

.551 

6.50 

.240 

.460 

.7 

.586 

5.57 

.182 

.500 

.8 

.625 

4.92 

.141 

.545 

.9 

.666 

4.45 

.110 

.590 

1.0 

.712 

4.10 

.090 

.620 

1.2 

.801 

3.60 

.065 

.670 

1.4 

.880 

3.23 

.055 

.690 

1.6 

.929 

2.98 

.049 

.710 

1.8 

.955 

2.80 

.047 

.715 

2.0 

.971 

2.68 

.046 

.720 
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TABLE n.- PARAMETERS IN DRAG COEFFICIENT EXPRESSION 


[From ref. 10] 


Mr 

Cd,C 

0.5 

0.48 

.6 

.455 

.7 

.497 

.8 

.535 

.9 

.610 

1.0 

.820 

1.1 

.900 

1.2 

.990 

1.3 

.993 

1.4 

.995 

1.5 

.999 

1.6 

1.000 

1.7 

1.020 

1.8 

1.005 

1.9 

.990 

2.0 

1.045 


Cd,fm 

A 

N 

7.80 

0.315 

0.410 

6.50 

.240 

.468 

5.57 

.182 

.500 

4.92 

.141 

.545 

4.45 

.110 

.590 

4.10 

.090 

.620 

3.85 

.070 

.645 

3.60 

.065 

.670 

3.42 

.060 

.680 

3.23 

.055 

.690 

3.11 

.052 

.700 

2.98 

.049 

.710 

2.89 

.048 

.713 

2.80 

.047 

.715 

2.74 

.047 

.718 

2.68 

.046 

.720 
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TABLE m.- PARAMETERS IN PRESENT DRAG COEFFICIENT EXPRESSION 


Mr 

Cd,c 

Cd,FM 

A 

N 

0.1 

0.380 

53.541 

1.7269 

0.1976 

.15 

.381 

35.759 

1.4099 

.2196 

.2 

.390 

26.888 

1.1908 

.2399 

.25 

.392 

21.580 

1.0339 

.2562 

.3 

.398 

18.053 

.9144 

.2706 

.35 

.403 

15.544 

.8159 

.2846 

.4 

.410 

13.670 

.7356 

.2973 

.45 

.419 

12.222 

.6672 

.3097 

.5 

.426 

11.065 

.6085 

.3215 

.55 

.435 

10.125 

.5637 

.3301 

.6 

.443 

9.345 

.5244 

.3384 

.65 

.453 

8.688 

.4890 

.3467 

.7 

.466 

8.128 

.4602 

.3536 

.75 

.480 

7.645 

.4367 

.3585 

.8 

.500 

7.224 

.4163 

.3630 

.85 

.513 

6.853 

.4043 

.3620 

.9 

.540 

6.525 

.3909 

.3631 

.95 

.600 

6.233 

.4735 

.3096 

1.0 

.710 

5.970 

.4384 

.3086 

1.1 

.780 

5.517 

.4332 

.3059 

1.2 

.820 

5.141 

.4261 

.3036 

1.3 

.860 

4.823 

.4252 

.3003 

1.4 

.890 

4.551 

.4260 

.2969 

1.5 

.910 

4.316 

.4334 

.2895 

1.6 

.920 

4.110 

.4392 

.2826 

1.7 

.930 

3.930 

.4483 

.2747 

1.8 

.940 

3.771 

.4535 

.2696 

1.9 

.940 

3.630 

.4545 

.2649 

2.0 

.940 

3.505 

.4489 

.2640 
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TABLE IV.- GAS FLOW -FIELD PARAMETERS FOR TEST CASES FOR VELOCITY -LAG CALCULATIONS 

SHOWN IN FIGURES 5 TO 7 


Type of flow 

PO 

T0» 

K 

Vg,x. 

m/s 

m/s 

Mr, max 

^Re,r,max 

Location of comparison 

lo5 N/m2 

atm 

Moo = 1.6, normal shock 

1.013 

1 

283 

218.4 

0 

0.7 

118 

1.27 cm behind shock 

Moo = 3, nozzle center line 

1.013 

1 

293 

616.8 

0 

.8 

35 

7.62 cm downstream of throat 

Moo = 3, normal shock 

1.013 

1 

293 

159.9 

0 

1.4 

10 

1.27 cm behind shock 

Moo = 3, 50 oblique shock 

1.013 

1 

293 

596.6 

52.5 

.3 

8 

1.27 cm behind shock 

Moo = 3, 10° oblique shock 

1.013 

1 

293 

564.1 

59.4 

.3 

12 

1.27 cm behind shock 

Moo = 5, nozzle center line 

3.446 

3.4 

363 

783.4 

0 

.9 

90 

54.61 cm downstream of throat 

' Moo = 5, normal shock 

3.446 

3.4 

i 363 

156.7 

0 

1.7 

1 56 

1 

1.27 cm behind shock 

Moo = 5, 50 oblique shock 

: 3.446 

; 3.4 

; 363 

■ 765.7 

; 66.7 

.4 

i 8 

1.27 cm behind shock 

Moo = 5, IQO oblique shock 

; 3.446 

3.4 

363 

^ 736.9 

129.9 

.7 

j 18 

i 1.27 cm behind shock 

j 

Moo = 6, nozzle center line 

3.446 

' 3.4 

, 554 

' 997.4 

0 

.8 

45 

j 168.6 cm downstream of throat 

Moo = 6, normal shock 

3.446 

' 3.4 

j 

554 

189.2 

0 

1.7 

! 

i 1.27 cm behind shock 

Moo = 6 , 50 oblique shock 

3.446 

1 3.4 

554 

977.4 

85.6 

.5 

' 3 

t 

1.27 cm behind shock 

1 

Moo = 6, IQO oblique shock 

1 3.446 

3.4 

j 554 

1 945.2 

166.7 

.9 

i 6 

1 

j 1.27 cm behind shock 


TABLE V.- CENTER-LINE NOZZLE DIMENSIONS 


Center -line nozzle flow 

Moo = 3 

8 

11 

CJl 

Moo = 6 

Distance from upstream injection point to throat, cm .... 

0.0 

15.24 

9.2 

Distance from throat to nozzle exit, cm 

7.62 

30.48 

168.6 i 

1 



Relative Mach number, 



O Aroesty (ref. 23) 

□ Sreekanth (ref. 24) 

O Ashkenas (ref. 25) 

A Zarin (ref. 26) 

Cii Coin and Lawrence (ref. 27) 
A Bailey and Hiatt (ref. 22) 


Figure 1.- Summary of available low Reynolds number and low Mach number experimental sphere drag data. 



Drag coefficient, 


O Bailey and Hiatt (ref. 22) 

□ Charters and Thomas (ref. 28) 

May and Witt (ref. 29) 
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(a) High Reynolds number (lO^ to 10®). 

Figure 2.- Comparison of available experimental sphere drag data. 
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(b) Mr = 2. 
Figure 2.- Continued. 
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(c) Mr = 0.2. 
Figure 2.- Concluded. 
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Figure 3.- Concluded, 
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Figure 4.- Comparison of present Cp) method with experimental data. 
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Cuddihy et al. (refs, 4 and 5) 

Carlson et al. (ref. 6) 

Crowe (ref. 7) 

Present method 



Dp, |j.m Dp, /im 


(c) Moo = 3; normal shock. (d) Moo = 3; 10° oblique shock x-component. 

Figure 5.- At the points noted in table IV, percent velocity lag predicted 

by various methods. 



3; IQO oblique shock y-component. (f) M«, = 3; 5° oblique shock x-component 



Dp, pm Dp, pm 

: 3; 50 oblique shock y-component, (h) Moo = 5; nozzle center line. 

Figure 5.- Continued. 




5; lOO oblique shock y-component. (1) Moo = 5; 5° oblique shock x -component 

Figure 5,- Continued. 




“ Cuddihy et al. (refs. 4 and 5) 

Carlson et aL(ref. 6) 


Crowe (ref. 7) 
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Dp, /im Dp, urn 

(m) Moo = 5; 5° oblique shock y-component. (n) Moo = 6; nozzle center line. 
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Dp, fJ-m Dp, fjLw. 

(o) Moo = 6; normal shock. (p) Moo = 6; oblique shock x-component. 

Figure 5.- Continued. 
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_ Cuddihy et al. (refs. 4 and 5) 

Carlson et al. (ref. 6) 

Crowe (ref. 7) 

Present method 




Moo = 6; IQO oblique shock y-component. (r) Moo = 6; 5° oblique shock x-component. 



Moo = 6; 5° oblique shock x-component. 


Figure 5.- Concluded. 



Crowe et al. (ref. 8) 

Waldman (ref. 9) 

Korkan et al. (ref. 10) 

Present method 



0 2468 10 0 2468 10 


Dp, ij.m Dp, iim 

(c) Moo = 3; normal shock. (d) Moo = 3; 10° oblique shock x -component. 

Figure 6.- At the points noted in table IV, percent velocity lag predicted 

by various methods. 
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Crowe et al. (ref. 8) 

Waldman (ref. 9) 

• • Korkan et al. (ref. 10) 

Present method 




(e) Moo = 3; 10° oblique shock y-component. (f) Moo = 3; 5° oblique shock x -component. 



Figure 6.- Continued. 
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Crowe et al. (ref. 8) 

_ . _ . _ Waldman (ref. 9) 

Korkan et al. (ref. 10) 

Present method 



(i) Moo = 5; normal shock. (j) Moo = 5; 10® oblique shock x -component. 



0 2468 10 0 2468 10 


Dp, pm Dp, pm 

(k) Moo = 5; 10® oblique shock y-component. (1) Moo = 5; 5® oblique shock x -component. 

Figure 6.- Continued. 
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Crowe et al. (ref. 8) 

- Waldman (ref. 9) 

• » * ' Korkan et al. (ref. 10) 
Present method 




Moo = 6; IQO oblique shock y~component. (r) Moo = 6; 5^ oblique shock x-component. 



Moo = 6; 5^ oblique shock y-component. 


Figure 6.- Concluded. 



Ill niHi ■■■1111111 iiiMi 


— Stokes Cj) equation (ref. 3) 

— Torobin and Gauvin (ref. 13) 

— Present method 
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(c) Moo = 3; normal shock. (d) M« = 3; 10° oblique shock x-component 

Figure 7.- At the points noted in table IV, percent velocity lag predicted 

by various methods. 



Stokes Cj) equation (ref. 3) 

— • — * Torobin and Gauvin (ref. 13) 

Present method 




3; IQO oblique shock y-component, (f) Moo = 3; 5° oblique shock x -component. 




3; 5° oblique shock y-component. (h) Moo = 5; nozzle center line. 


Figure 7.- Continued. 




Stokes Cj) equation (ref. 3) 
Torobin and Gauvin (ref. 13) 
Present method 




(m) Moo = 5; 5° oblique shock y-component. (n) Moo = 6; nozzle center line. 




Moo - 6; 10*^ oblique shock x -component. 


Dp, /Liin 

(o) Moo = 6; normal shock. (p) 

Figure 7.- Continued. 
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(a) Particle size, 1 jum. 



(b) Particle size, 3 jum. 



Stokes Cd equation (ref. 3) 
Present method 


(c) Particle size, 5 jUm. 

Figure 8.- Variation of percent velocity lag behind a Mach 1.6 normal shock (stagnation 
conditions listed in table IV) predicted by various methods. 
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(a) Particle size, 1 Mm. 



(b) Particle size, 3 /xm. 



Stokes Cj) equation (ref. 3) 
Present method 


(c) Particle size, 5 Mm. 

Figure 9. - Variation of percent velocity lag behind a Mach 3 normal shock (stagnation 
conditions listed in table IV) predicted by two methods. 
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(a) Particle size, 1 jum. 



X2, cm 


(b) Particle size, 3 /im. 



Stokes Cj) equation (ref. 3) 
Present method 


(c) Particle size, 5 /xm. 

Figure 10.- Variation of percent velocity lag behind a Mach 5 normal shock (stagnation 
conditions listed in table IV) predicted by two methods. 


53 



(a) Particle size, 1 jum. 



(b) Particle size, 3 fxm. 



X 2 , cm 


Stokes Cd equation (ref. 3) 
Present method 


(c) Particle size, 5 ixm. 

Figure 11.- Variation of percent velocity lag behind a Mach 6 normal shock (stagnation 
conditions listed in table IV predicted by two methods. 
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